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Abstract 

A recently proposed scheme for the analytical treatment of the dynamics of 
two-dimensional hybrid Josephson junction arrays is extended to a class of 
generalized hybrid arrays with "horizontal" shunts involving a capacitive as 
well as an inductive component. This class of arrays is of special interest, 
because the internal cell coupling has been shown numerically to favor in- 
phase synchronization for certain parameter values. As a result, we derive 
limits on the circuit design parameters for realizing this state. In addition, 
we obtain formulas for the flux-dependent frequency including flux-induced 
switching processes between the in-phase and anti-phase oscillation regime. 
The treatment covers unloaded arrays as well as arrays shunted via an external 
load. 
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I. INTRODUCTION 



Two-dimensional Josephson junction arrays are considered as strong candidates for tun- 
able microwave oscillators. Since the pioneering works by Benz and Burroug hHi there 
were some attempts to fabricate arrays of this typeH as well as to understand them 
theoreticallyiHi. (For two recent reviews on 2D Josephson junction arrays, see Lachenmann0 
and Booi0.) While radiation output of two-dimensional arrays should be much larger than 
that of one-dimensional arrays (for quadratic arrays in the matched case typically ~ iV 2 , 
with iV=number of rows, compared to being ~ N, with iV=number of junctions, for one- 
dimensional arraysi'0) observations point more to the opposite direction. While in one- 
dimensional arrays there were observed up to 160/iW, 

0.0 

the output power reported in 

two-dimensional arrays is several orders of magnitude smaller with a maximum of around 
400nW§S§. 

Potentially, there can be several reasons responsible for this discrepancy. Besides low 
critical currents/normal resistances, mismatch to the external load, or parameter tolerances 
some more basic problems might be responsible for this. Indeed, some recent theoretical 
investigations show, that the radiating in-phase mode is neutrally stable in a unshunted array 
without external fluJ and, even worse, that it is unstable even for a small flux entering the 
celli. As a result, the natural state of at least the simple model circuit studied inEH is 
a non-radiating one with both cells oscillating against each other. The situation can be 
improved by adding an appropriate external shunt synchronizing in-phase via its long-range 
interaction, but generally there remains a tendency that pairs of cells lock anti-phase and 
drop out of the radiating mode. 

A recently proposed layout0 removes this difficulty by introducing an additional ca- 
pacitive shunt in the "horizontal" branches thus turning the internal coupling to favor the 
in-phase state. It is the aim of the present investigation, to give this idea an analytical 
foundation and derive some rigorous results, notably on the parameter boundaries sepa- 
rating in-phase from anti-phase oscillations. In addition, we study the interplay with an 
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external load leading to a rather complex picture of possible stability regions as a result of 
the competition of external and internal synchronization. 

We start with an exposition of the problem including the basic equations in Sec. II. In 
Sec. Ill and IV these equations are solved by an analytical approximation, combining ideas 
of the strong coupling method appropriate for small-inductance Josephson junction cells0 
with the standard weak-coupling procedure of slowly- varying phaseliU^ for the treatment of 
inter-cell coupling. While Sec. Ill contains lowest-order results corresponding to vanishing 
cell inductance, Sec. IV includes the effects caused by a small, but non- vanishing inductance 
being essential for understanding the inter-cell coupling. Sec. V contains several results 
including a comparison with numerical simulations. The interplay with an external load is 
treated in Sec. VI, and Sec. VII contains several more general conclusions relevant for the 
layout of two-dimensional Josephson junction arrays. 



II. THE MODEL AND THE BASIC EQUATIONS 

For making the problem accessible to an analytical treatment, we have to make several 
propositions. Fig. 1 shows the circuit under consideration. To make the physical mechanisms 
more transparent, the external shunt Z$ is removed in the beginning, and will only be 
included in Sec. VI. Despite its simplicity this model has all the main features present in 
larger arrays, too: It is truly two-dimensional with a possible flux entering the cells. Notice, 
that unlike to conventional hybrid arraysii the horizontal branch contains a more general 
shunt consisting not only of the usual inductive connection, but of a parallel capacitance 
and resistance. Numerical results obtained before indicate that an intrinsic shunt of this 
type can favor in-phase locking even in an externally unshunted array. Here, we will confirm 
and extend this result by developing an analytical formalism which should be applicable to 
larger two-dimensional arrays with the same general structure as well. 

Some more restrictions have to be put on the array: (i) Josephson junctions are described 
by the RSJ mo (ii) All junctions are considered to be identical, (iii) Junctions are 
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overdamped with a McCumber parameter set to zero, (iv) Self-inductance is taken into 
account while mutual inductance is neglected, (v) The normalized ring inductance between 
the two loops 

I = 2nI c L/% (1) 

is supposed to be small (7 <C 1). From the beginning, one has to understand that the 
inductance of the horizontal connection acts in a two-fold way. At first, it contributes to 
both ring inductances thus determining the SQUID coupling within each loop. At second, it 
is part of the shunt common to both loops and as such it influences the inter-cell coupling. 
With Eq. (P we request the SQUID coupling to be strong, which is a necessary prerequisite 
for our approximation scheme to work. On the other hand, we will not fix the ratio between 
the inductive and the capacitive horizontal impedances from the beginning. 
In the following, we will exploit some more normalized quantities as follows, 



s = —I c R N t, (2) 

y? = 27T$/$ , (3) 

c = ^I C R 2 N C, (4) 

r = R/R N , (5) 

i = I lie, (6) 



with Ic the junction critical current, R^ the normal resistance of one of the (identical) 
junctions, $ the external flux per cell, and the last normalization being valid for all currents 
entering the calculation. Adopting these normalizations, the circuit can be described by the 
RSJ equations for the Josephson phases, 

fa + sin0y = it,, ({i, j} = {1,2}) (7) 

in conjunction with the two flux quantization conditions 

<fri2 -(j>a-cp^lii = (8) 
5 



(minus sign refers to % = 1) and Kirchhoff's current laws 



h = ^(iii + 112), (9) 
« = in - hi = 122 - «i2, (10) 
ii = i + i rc . (11) 

These have to be supplemented by Kirchhoff's voltage law 

T • I 

ii + -jhc + -j^hc = 0. (12) 

We would like to point out that while the inductive branch carrying current i\ is part of 
both superconducting loops thus contributing to the flux quantization conditions the 
branch i rc enters only via the ordinary Kirchhoff's law flT2]). As a result, it is impossible to 
simply substitute the three elements /, c and r by a single impedance Z from the beginning. 

Before, it has proven useful in the treatment of strongly coupled SQUID ce to combine 
the Josephson phases within each cell via 

s fc = - {4> k 2 + 4>ki) , (13) 

A fc = i {(t) k2 - 0fci) • (14) 

In addition, we introduce the circular currents 

ifc = (i*2-i fc i)/2. (15) 

With the help of Eqs. ( |13|) - dl5|) we finally obtain the system 

Sfc + sin Efc cos = i , (16a) 

A fe + sin A fc cosS fc = i° k , (16b) 

A 1 + A 2 -<^ = 0, (16c) 

Ai — A 2 = Hi — l{%2 — i rc ), (16d) 



r • 1 

Vc + tVc + -r^rc = (h ~ h°) (16e) 
I Ic 

which our analytical approximation scheme is based on. As there are seven equations for 
the seven variables E fc , A k , i° k) i rc , this is a well-posed problem. 
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III. ANALYTICAL APPROXIMATION SCHEME AND LOWEST ORDER 

RESULTS 



Our strategy for solving system ( |T5] ) will be based on a perturbative treatment valid for 
small I (for the basic idea compare our earlier paperEl). Thus, we start solving Eqs. ( |To] ) 
for Z = 0, and only later include corrections ~ / exploiting the lowest order results obtained 
before. This procedure is favored by the fact that I enters equation (|16d|) only We start by 



evaluating (|16d|) in conjunction with (|16c|) . The solutions for A& can be used to evaluate 
Efc from ( |16a| ). Next, we find the i° k (not the A*., which are already known in this order!) 
from (|16b| ). Finally, with the ring currents i° k on the right hand side of (|16e|) known we 



can evaluate the current i rc by solving the corresponding differential equation. All other 
quantities, like %i or i, are secondary and can be derived from the seven variables mentioned 
so far. Afterwards, we insert the lowest order result on the right hand side of Eq. (|16d| ) and 
start a second cycle in the same sequence. 

The procedure described above gives the following lowest-order results. First, the Joseph- 
son phase differences in both loops are found to be identical, 

A M = <p/2. (17) 



In the following, comma-delimited indices refer to the order of approximation. From 
the Josephson phases are found to coincide with the corresponding solutions for an au- 
tonomous junction,^! 

£ fc ,o = 77 + 2arctan ( — - — — tan ^ 5k ) , (18) 
2 \i + cos(y?/2) 2 J 

with the important modification that the frequency Co becomes flux-dependent according to 

Co = ^§-cos 2 (p/2). (19) 

Next, the circular currents can be evaluated from ( |16b| ). Note, that this equation origi- 
nating from the original Josephson equations does not lead to a differential equation, because 
the constant Josephson phase differences A^ are already known. The result is 



i° k = sin(y?/2) cosE M . (20) 



It is a trivial task to evaluate Eq. ( |20"D using (fHf); in the further calculation we will only 
need the lowest harmonics of the circular currents, 

il,o = sin (^ 2 ) sin (Co^ - S h ). (21) 

^0 + C,0 

The corresponding difference of the ring currents, 

enters the horizontal connection thus acting as a driving force for the oscillatory circuit 
according to Eq. (|16ej ). This equation can be solved with standard methods. The stationary 
oscillating solution reads 



4/Co 2 , , /nN • (Si-6 2 \ . , 5i + 5: 



irc '° = "|z(Co)Kio + Co) sin(¥,/2) sin l^V 1 J sin l CoS " ^ - ^ (Co) J • (23) 

Here, we introduced the series circuit impedance Z with 

\Z(Co)\ 



r 3 + ( 4r - ^Co ) (24) 



and the phase angle -0, 



"■^ = izfer- 8in ^» )= ikcf- (25) 

For later purposes we need ^ = i+z rc rather than z rc , because it is just i\ which potentially 
may split the oscillation phases between cell 1 and cell 2 via Eq. (|16d ). Combining fl23|) 
with (|22|) after some algebra we obtain 

...o = sm cos (Co* - - x(Co) j , (26) 

where we introduced the rc impedance \z\ with 



W6)l = f 2 + w (27) 
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and 



, r 2 r " c (o +(^r- Ko) 

smx(Co) = ^ cos X (Co) = P° ' • (28) 

Z(Co I Jl + (rcCo) 2 \Z(Co) | Jl + (r<o) 2 



(In principle, one could evaluate i\ directly from an equation similar to ( |16e| ) of course, and 
we checked that the result is the same. The procedure described here has the advantage of 
additionally providing an expression for the current flowing through the capacitive line.) 
To summarize, we observe the following lowest-order results: All four junctions oscillate 



with the same flux-dependent frequency Co = \Ao — cos 2 (<^/2). Because of fll7D , the junctions 
within each cell are exactly in phase, while the relative phase between cell 1 and cell 2 (given 
by Si and <5 2 , respectively) is undetermined, up to now. If both cells are in phase, there is 
no current through the horizontal line, because of the sin[(<5i — 62) /2] present in (|22|) . On 
the other hand, the horizontal current reaches its maximum if both cells oscillate anti-phase 
with Si — 62 = 7r. 

IV. INDUCTANCE EFFECTS 



Now we are ready to include inductance effects. Again starting with ( |16c| ) and ( 16d ), we 
insert the lowest-order result ( p^) on the right hand side of (|16d| ). This leads to 

Ai /2 = Ai /2 ,o + lAi /2A (29) 
V , , Co \z\ . , , , . S 2 -Si f S1 + S2 \ 



Note, that the first index in ( P9|) refers to cell 1 and cell 2, respectively, while the second 
one indicates the order of evaluation; the + sign refers to A x . This has to be inserted into 

(Hi), 

t k + cos(A M + lAk,i) sin S fc = i . (30) 
For evaluating these equations the cosine on the left hand side is expanded according to 

cos( A fcj0 + lA k ,i) ~ cos A M - lA k ,i sin A fc)0 . (31) 



After transferring the correction term ~ / to the right hand side of (pOj) one makes the crucial 
observation, that it acts in a similar way as,for example, an external shunt synchronizing 
the cells0. 

The resulting equations are evaluated with the conventional phase-slip method (see, for 
instance&~§3). According to this procedure which has proven useful in the treatment of linear 
arrays before, the up to now constant phases 5\ and 62 are considered as time dependent, 

4 = S k (s), (32) 

with the subsidiary condition that this time-dependence is only an adiabatic one, 

5 « Co- (33) 

Physically, this means that the phases are required to be nearly constant during one Joseph- 
son oscillation. 

With these assumptions, the same ansatz flUf ) with 5k{s) and C instead of Co leads to the 
sum voltages 

£ k = Co(C - h) ^ 

i + cos(</?/2) cos(Cs — 5k) 

Writing ( instead of Co we have allowed for a possible (small) deviation of the actual oscil- 
lation frequency from Co- Inserting fl3l| ) into (|30| ) leads to the reduced equations 

Co(C - Co - h) = Zsin(<^/2)A M (cos(y9/2) +i cos((s - S k )) . (35) 

After averaging over one time period and applying some algebra we arrive at the following 
system of equations (for details see, for instance@'&il) 

Co(C - Co- < $1 >) = lio . t4t sin(y?/2) sin < ^ > < - > sin(<p/2) 

ia + Co\Z\ 2 

(< 5 2 > - < 81 > \ ( . 

x cos I h X I , (36a) 

Co(C - Co- <S 2 >) = -li ■ ^° , t|ft sin(v9/2) sin < ^ 2 > < - ~ h > sin(y?/2) 

*o + so 1^ I 1 

(<5 2 >-<6 1 > \ 
x cos h x , (36b) 
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where <> denotes the time average over one Josephson oscillation. The difference of (|36aj) 
and ( |36b| ) gives an evolution equation for the phase difference < 5 >, 



< 5 >= . %q1 - i ^ i -sin 2 (^/2) cos x sin < 5 > . (37) 
z o + Co \Z\ 

Eq. ( P?D is the basic equation determining the possible phase differences between the oscil- 
lations of both cells as well as the corresponding regions of stability. 

V. PHASE LOCKING, STABILITY AND OSCILLATION FREQUENCY 

We will not go into the question of general solutions of ( |37|) but concentrate on phase- 
locking, being characterized by a time-independent phase-shift between cell 1 and cell 2, 

< S lock >= 0. (38) 



Within the range < 5 < 2tt there are obviously only two possibilities for Eq. ( p8|) to be 
valid, 

<5 lock >=0 and <5 lock >=7r, (39) 

the first one describing in-phase oscillations and the second one anti-phase oscillations of 
the cells. 

The crucial question of the range of stability of these two solutions can be answered on 
the basis of the evolution equation fl3~7|), too. The ansatz 

< 5 >=< 5 lock > +ae xt (40) 

(\a\ |1|) leads to the Lyapunov coefficient 

A = -i^_M sin 2(^/2) cosxcos < s lock > . (41) 
l o + Co \Z\ 

One recovers, that the stability is solely determined by the cos x] all the remaining factors, 
except 5 lock , are positive definite. In detail, the 
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in-phase solution < <5 lock >= is stable for cosx < 0, (42) 

while the 

anti-phase solution < <5 lock >= tt is stable for cosx > 0. (43) 

Before further evaluating this condition we will consider the oscillation frequency which 
can be derived from ( |36a| ) (or ( |36b|) ). With 



< St >=< 5 2 >= const = (44) 



one easily recovers 



C in = Co = v^-cos 2 ^). (45) 
Evaluating the anti-phase frequency with 

< 5i > - < 5 2 >= n (46) 
needs a bit more algebra. The result is 



Thus, if both cells oscillate in-phase their frequency is identical to the autonomous oscillation 
frequency. On the other hand, if the cells oscillate anti-phase the frequency will be lower 
than Co- The physical reason for this behavior can be understood by comparing with other 
(even linearly) oscillating systems: If the bindings (in our case realized by the horizontal 
impedance) are not loaded, the oscillation frequency remains the same as for uncoupled 
oscillators; if the bindings are loaded (i.e. in case of a ac current flowing through the 
horizontal line) the system oscillates with a different frequency. 



Unfortunately, one has to respect a certain limit of validity of Eq. fl4"7). Using the 
method of slowly varying phase we have adopted the supposition mentioned before that the 
frequency must not deviate too much from Co, 

C ~ Co- (48) 
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Thus, the correction in Eq. ( [47]) is required to be small compared to the frequency itself. A 
rough estimate valid for i > 1.15 leads to the condition 



Our experience shows that usually a factor of 2 ... 3 is sufficient for this condition to be 
fulfilled. 



Now we return to the question of anti-phase «-> in-phase transitions described by (|42|) and 
([43|), resp. Considering the numerator of cosx one observes that the boundary separating 
in-phase and anti-phase oscillations of the cells is given by 



with the cells oscillating anti-phase if the left hand side is positive and in-phase if it is 
negative. In other words, the transition between both regimes lies in the vicinity of the 
resonance curve of the l-c-r connection with deviations becoming important for small I. Fig. 
2 shows the boundary between the two regimes for a frequency ( = 1.11 in comparison to 
numerical results. 

To summarize, the in-phase regime is favored for not too large r as long as the inductive 
impedance dominates over the capacitive one, while for the capacitive impedance dominating 
the cells oscillate anti-phase. There is a simple physical explanation for this: Anti-phase 
oscillations are caused by the flux coupling via the joint inductive line carrying current i\. 
For a sufficiently large capacitive shunt, the current prefers the capacitive way which does 
not produce any such flux. 

In conventional hybrid arrays! horizontal lines are purely inductive. Formally, this limit 
can be observed letting c — > 0. In this case the capacitive impedance goes to infinity while 
the correction ~ r 2 c tends to zero. Then, there is no possibility for the current to be 
shunted, and the cells remain in the anti-phase regime@. The more general question, for 
which parameter values I, c, and r there are no transitions can be answered on the basis of 
Eq. d50|) . This equation does only have real solutions for ( if 



(49) 




(50) 
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/ > r 2 c. 



(51) 



For all smaller I, the current in the inductive line is strong enough to keep the cells oscillating 
anti-phase. 

Considering the circuit parameters io, I etc. as constant and leaving the external flux (f 
as the only free parameter one can observe flux-induced transitions between both regimes. 
The difference between the frequencies ( m and ( antl leads to a hysteresis, which has been 
observed in numerical simulations beforeS. In more detail, in-phase — > anti-phase transitions 
are observed at 



<P 



'' — 2 arccos 



tr2 



(52) 



where we introduced the transition frequency 



c 



tr 



1 



ylc — r 2 c 2 



(53) 



as can be easily deduced from flSOl). The transition from the anti-phase to the in-phase 
regime needs a bit more algebra. It can be determined from the requirement, that the anti- 
phase frequency ( fi7| ) be equal to the transition frequency (|53"D. Unfortunately, the resulting 
equation can not be solved in closed form. However, as a first approximation, one can equate 
the in— >anti transition frequency ( (5"3"D with and evaluate for if, substituting ( — > ( tr on 
the right hand side of (|47|) , 

cr U + C 

± kA • (54) 



2 arccos 



V 



\ 1 cr + C tr ) 



It can be deduced, that <p ai is always larger than ip m . A better result for <p ai is obtained by 
graphically finding the transition frequency on the curve at C — C tr - 

Thus, if there are any transitions between both regimes at all, for small values of the 
external flux the cells oscillate with the lower anti-phase frequency switching to in-phase 
oscillations at y? ai . Because of ip ai > ip ia (for < if < n) switching back to the anti-phase 
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state occurs at a lower flux, leading to the hysteresis mentioned above. Fig. 3 shows a plot 
of frequency against flux in comparison with the outcome of a numerical simulation. The 
frequencies are in excellent agreement, and even the transition points, which depend rather 
sensibly on the parameters, are located within the same region. 

This last result concerning hysteresis has to be taken with some care. It was obtained 
by combining the anti-phase frequency formula fl47p with Eq. (|50|) and evaluating for (p. 



However, (|50| ) as originating from fl37|) is already a first order result, thus inserting 
might not be fully justified while second order terms in ( |37l) are neglected. Nonetheless, it 
gives a plausible explanation for the mechanism causing the hysteresis observed in numerical 
simulations. 

VI. LONG-RANGE SYNCHRONIZATION VIA AN EXTERNAL LOAD 

It has been well-known for a long time that synchronization in a one-dimensional array 
can be achieved and controlled by shunting the array via an external load0@. In a similar 
manner one may hope to be able to control row locking in two-dimensional arrays, too. For 
studying this mechanism within our model we now add the external load already indicated 



in Fig. 1. As a result, we have to supplement the basic equations ([Tfj). At first, we add the 
mesh rule for the load current i L , 

S fc - hh ~ r L k - —ih = 0. (55) 

k=l,2 ° L 

Here, vl,Il, an d cl are the load impedances normalized in the same manner as ([!]), and 
(§). In addition, the load current couples back to the junctions, thus supplementing Eq. 

E fc + sin£fcCos A k = i - ^i L . (56) 

As has been observed in the study of similar one- dimensional synchronization problems 
before, the reciprocal impedance \ j\Zj\ <C 1 provides another perturbation parameter for a 
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sufficiently large load; thus we evaluate the system perturbatively, neglecting terms ~ l/\Zi\. 
To lowest order with respect to \Zj\ the load current vanishes, and we end up with the 
results described in Sec. III. Based on the lowest order Josephson oscillations (|I5D and the 
corresponding voltages E^o we obtain the first order (with respect to 1/\Zl\) load current 

iL,o, 

4cos(y/2) Co f Sx-SA . / 5 1 + 5 2 , \ 
il,o = n=-j ■ , , cos — - — sin Cos ip L (57) 



io + Co 



with 



l^(Co)| 



(r L + I) 2 + ( —7- - kCo] , (58a) 



V L ' ri ' "V^Co 



sin ip L (Co) = ,y ^tt , (58b) 
' /lCo 



cos^(Co)= Cj *° ^, • (58c) 

Its structure is obviously quite similar to that of the horizontal current (|22|). However, one 
should note two differences: (i) While the load current is maximal for tp = 0, the horizontal 
current reaches its maximum for (p = tt/2. (ii) The horizontal current vanishes if both cells 
oscillate in-phase, while the load current vanishes for both cells oscillating anti-phase. 

The load current fl57D provides the additional contribution to (|56|) and, as a result, the 
phase slip equations (^) get an additional term, too. After performing the time-averages 
we get 



(59a) 



Co(C - Co- <6 1 >) = -4^TTW\ sin2 (¥V 2 ) + sin(< 5, > - < 5 2 > - X )) 

2 i + Co \Z\ 
1 C 1 

+ - . ° 1 „ 1 cos 2 (yV2) (sin^ £ - sin(< S 1 > - < 5 2 > -ip L )) , 
2i + Co \Z L \ 

Co(C - Co- <S 2 >) = ^-^-M s i n ^/2) (sin ^ + sin(< 5 1 >-<6 2 > - X )) 

+ 77 • ,°> , 7 1 cos 2 (v?/2) (sin^ L + sin(< ^ > - < 5 2 > -^l))- 
^ ^o + c,o 

By subtracting ( |59a| ) and Q59b| ), we finally get the evolution equation for the averaged 
oscillation phase difference, 
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< S >= -. — (r^—r cos 2 (ip/2) cosipL + ^o^rjf sin 2 (y?/2) cosx J sin < 5 > . (60) 

io + Co \\Zl\ \Z\ J 

Despite the relatively complicated interplay between cell interaction via the horizontal line 
and long range coupling via the external load there remain only the same two phase locking 
solutions as before, 

<5 lock >=0 and <5 lock >=7r, (61) 

the stability of which is determined by the Lyapunov coefficient 

A = — - — ( cos 2 (^2) cos ^ £ + i J M sin 2 Up 12) cos X ) cos < <5 lock > . (62) 
i + Co \ \Zl\ \Z\ J 

In-phase oscillations of the cells are stable if the term in parenthesis is lower than zero 
while anti-phase oscillation are stable if it is greater than zero. Thus, the desired stability 
for the in-phase mode is reached for 

1 \z\ 

■ cos 2 ((^/2) cosipL + ^7777 sin 2 ((£>/2) cosx < 0. (63) 



\Z L \ " vr/ ' r " \Z\ 

Eq. ( pB]) shows a rather complex parameter dependence, relating the seven parameters 
r,l,c,TL,lL,CL, and <p. Its physical meaning is best discovered considering several limiting 
cases. 

(i) For a sufficiently large external load, 

1 ( 64 ) 



\Zl\ \Z\ 

the relative phase of the cells is determined by the internal coupling alone. This has to be 
compared to the case of two externally loaded separate cellsil. In this CclSG clS for linear 
arrays - the relative phase depends on the character of the external load only: While for 
inductively dominated loads the cells are locked in-phase, they are locked anti-phase for 
capacitively dominated loads, independently of the magnitude of the external load. 

(ii) The contributions from the external load and from the internal shunt show a different 
flux dependence. For sufficiently small values of external flux the last term can be neglected, 
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and the locking regime is controlled by the load only. On the other hand, for flux values of 
around half a flux quantum the first term becomes negligible, and the internal horizontal 
line determines the phase difference of the cells. 

(iii) For I — > 0, the second term can be neglected, and the result agrees with that obtained 
for two separate cells before^, as it should be. In this limit the cells internally decouple, 
while the external coupling remains in force. 

(iv) The usual hybrid arrays without the internal i?-C-line are contained as a limiting 
case. For r — > oo, the in-phase- condition, Eq. (|63|), reduces to 

—J— cos^ L cos(> 2 /2) + i l sin 2 0/2) < 0. (65) 

It states, that for sufficiently large inductances, 

l>l = ~ ■ I 2? ,n\\v v 66 

ordinary pure inductive hybrid arrays may switch to the anti-phase state even for inductive 
external loads. 

The indicated transition was indeed observed in a numerical simulation (boxes in Fig. 
4). Having in mind that Eq. (pq ) is the result of several approximations, concerning the 
external shunt as well as the internal inductive coupling, the agreement is remarkably good. 

The influence of changing parameters can be nicely illustrated by performing a second 
simulation with exactly the same parameter set, but distributing ring inductance I regularly 
around the loops. The result denoted by the dots in Fig. 4 clearly deviates from that 
obtained for inductance concentrated on the horizontal line considered before. This can be 
taken as a strong indication that the coupling is not provided by the loop inductances but 
by the inductance on the line common to both cells. 



VII. CONCLUSIONS 



Although our work is devoted to the study of a simple model circuit several results 
are expected to be valid for larger arrays, too. At first, the short range coupling between 
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neighboring cells leads to an anti-phase synchronization in conventional Josephson junction 
hybrid arrays. This may be one reason for explaining the very small radiation output in 
2D Josephson junction arrays obtained so far. At second, we show a way to improve the 
situation by adding a capacitive shunt parallel to the horizontal lines. In this way, the flux- 
generating current potentially responsible for the anti-phase coupling is redirected through 
the capacitive line which is not part of a flux quantization condition. 

Combining Fig. 2 with some already known facts on synchronization in strongly coupled 
SQUID cellsS the following design criteria for generalized hybrid 2D Josephson junction 
arrays can be derived, (i) For synchronizing horizontal lines in-phase the ring inductances 
have to be kept small (/ <C 1). (ii) In-phase synchronization between neighboring cells in 
vertical direction is observed for I > l/c( 2 + r 2 c. Based on this, we will derive some estimates 
for reasonable c and I. For a given I, the boundary between in- and anti-phase oscillations is 
given by Eq. (|50|). Fig. 2 shows already that the additional term ~ r restricts the possible 
I by setting a lower bound. This bound is obtained from 

as 

r = — resp. / = 2r/(. (68) 
Thus, for obtaining in-phase oscillations the condition 

r < | (69) 



has to be respected. Because of (|68|) this means 



c>£. (70) 

Obviously, the requirement to have a small I for horizontal in-phase synchronization leads 
to the demand to have a sufficiently high capacitance c > 2/7£ as well as a small resistance 
r < l(/2. A reasonable compromise might for instance be 
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Z«0.8 c^3.0 r^0.2. (71) 

Of cause, all these estimates should be considered as very rough, and on the other hand, 
one has to check carefully how large these quantities on chip actually are. 

On the other hand, we would like to point out that these suggestions are based on an 
analytical approximation scheme and are founded on solid formulae. Of cause, it still has to 
be shown rigorously that they can be transferred to larger arrays as well. Some preliminary 
results from numerical simulations indeed indicate this. We hope, that the general procedure 
described here can be transferred to larger arrays of the type considered here as well, and 
some work is on the way to actually extend it to a ladder configuration. 

If the arrays are externally loaded, which is usually done via an inductive load, the 
parameters have to be chosen in such a manner to respect Eq. fl63l). The best way for 
obtaining in-phase synchronization is to make both contributions to the Lyapunov-coefficient 
lower than zero separately, which is possible because the parameters of the external load 
can be chosen independently of those from the internal shunt. In general, one should select 
values such, that (i) the external load is dominated by its inductive contribution, (ii) the 
internal horizontal shunts are dominated by the inductive impedance, too. Because of the 
frequency-dependence of the characters of the shunts, one has to make sure, that these 
conditions are met for all values of external flux. 

Of cause, the circuit studied here has several features requiring a more detailed investi- 
gation, either analytically or numerically. Usually one exploits shunted tunnel junctions for 
building arrays, thus one may ask for the influence of non-vanishing McCumber parameters. 
On the other hand, the influence of parameter splitting needs to be investigated, and in 
addition, in real arrays, noise comes into play. While this last aspect is to be expected to 
play only a minor role within the small inductance loops, it will be sure have some influence 
on the coupling between the cells. 
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FIGURES 

FIG. 1. The generalized Josephson junction hybrid array model circuit under investiga- 
tion. 

FIG. 2. The boundary between in-phase and anti-phase oscillations. Solid line: analytical 
approximation. Crosses: numerical simulation. Parameters: io = 1.5, r = 0.1, <p = 1.0. 

FIG. 3. Frequency against flux with a transition from anti-phase to in-phase oscillations. Pa- 
rameters: io = 1.5, r = 0.1,1 = 0.2, c = 4.0. (a) Analytical approximation, (b) Numerical simula- 
tion. 

FIG. 4. Transition between in-phase and anti-phase state caused by the internal inductive 
coupling present in a hybrid array. Parameters: io = 1.5, tl = 1.0, Zt, = 1.0, cl = 2.0. Solid 
line: analytical approximation, boxes: numerical simulation, open dots: numerical simulation with 
inductance regularly distributed around loops. 



22 



REFERENCES 

1 S. P. Benz and C. J. Burroughs, Supercond. Sci. Technol. 4, 561 (1991). 
2 S. P. Benz and C. J. Burroughs, Appl. Phys. Lett. 58, 2162 (1991). 
3 P. A. A. Booi et al, IEEE Trans. Appl. Supercond. 3, 2493 (1993). 
4 L. L. Sohn et al, Phys. Rev. B 47, 975 (1993). 

5 H. Shea, M. Itzler, and M. Tinkham, Phys. Rev. B 51, 12690 (1995). 

6 K. Wiesenfeld, S. Benz, and P. Booi, J. Appl. Phys. 76, 3835 (1994). 

7 G. Filatrella and K. Wiesenfeld, J. Appl. Phys. 78, 1878 (1995). 

8 R. L. Kautz, IEEE Trans. Appl. Supercond. 5, 2702 (1995). 

9 M. Darula, P. Seidel, F. Busse, and S. Benacka, J. Appl. Phys. 74, 2674 (1993). 

S. Lachenmann, Ph.D. thesis, Eberhard-Karls-Universitat Tubingen, 1995. 

1 P. A. Booi, Ph.D. thesis, Twente University, 1995. 

2 K. K. Likharev, Dynamics of Josephson junctions and circuits (Gordon and Breach, 
Philadelphia, 1991). 

3 J. E. Lukens, in Superconducting Devices, edited by S. T. Ruggiero and D. A. Rudman 
(Academic Press, New York, 1990), pp. 135-167. 

4 J. A. Stern, H. G. LeDuc, and J. Zmudzinas, Trans. Appl. Supercond. 3, 2485 (1993). 
5 M. Octavio, C. B. Whan, and C. J. Lobb, Appl. Phys. Lett. 60, 766 (1992). 

6 M. Basler, W. Krech, and K. Platov, subm. to Journ. Appl. Phys. (1996). 

7 W. Krech and K. Platov, subm. to Europhys. Lett. (1996). 

8 M. Basler, W. Krech, and K. Y. Platov, Phys. Rev. B 52, 7504 (1995). 

23 



A. K. Jain, K. K. Likharev, J. E. Lukens, and J. E. Sauvageau, Phys. Rep. 109, 310 
(1984). 

W. Krech, Ann. Phys. (Leipzig) 39, 117 (1982). 
W. Krech, Ann. Phys. (Leipzig) 39, 349 (1982). 
W. Krech, Ann. Phys (Leipzig) 39, 50 (1982). 

M. Basler, W. Krech, and K. Y. Platov, in Macroscopic Quantum Phenomena and Co- 
herenece in Superconducting Networks, edited by C. Giovanella and M. Tinkham (World 
Scientific, Singapore, New Jersey, London, Hong Kong, 1995), pp. 225-234. 



24 




■\ -\ I I I I I I L 

0.5 1 1.5 2 2.5 3 

Flux (p 




■\ -\ I I I I I I L 

0.5 1 1.5 2 2.5 3 

Flux (p 



